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COMPUTER EN WISKUNDE !) 
door 


Prof. Dr. A. VAN DER SLUIS 
Utrecht 


1. De komst van de computer heeft de beoefening van verschei- 
dene takken van wetenschap op revolutionaire wijze beïnvloed. 
Ten aanzien van de zuivere wiskunde is dit nauwelijks het geval. 
Wel heeft men de computer daarbij o.a. gebruikt om groepen en 
algebra’s van zekere structuur te vinden; om groepkarakters uit 
te rekenen; om grote aantallen 10 x 10 grieks-latijnse vierkanten 
te construeren (Euler dacht dat er niet een bestond); en om bij 
het opkomen van een ad-hoc vermoeden snel een aantal gevallen 
door te rekenen: alvorens moeizaam naar een bewijs te gaan zoeken. 
Maar ook de nieuwste ontwikkelingen in de zuivere wiskunde zijn 
zonder toedoen van de computer tot stand gekomen, en de beroemde 
vermoedens (Fermat, Riemann, Goldberg, vierkleurenpro- 
bleem) bestaan nog steeds. 

2, In plaats van wiskundige problemen op te lossen heeft de 
computer veeleer nieuwe gesteld, en zo aanleiding gegeven tot de 
ontwikkeling van nieuwe mathematische disciplines en de opleving 
van oude. Voorts heeft de computer-gebruiker (vaak door schade 
en schande) een eigen visie ontwikkeld ten aanzien van de, vaak 
door de traditie geheiligde, voorraad conventionele wiskundige hulp- 
middelen. In verband met het thema van de dag (computer en 
onderwijs) zal met name dit laatste onze aandacht hebben. 

3. Een schoolvoorbeeld van een dergelijk hulpmiddel is de zgn. 
abc-formule voor de oplossing van een vierkantsvergelijking, een 
hulpmiddel waarmee ieder zich volkomen vertrouwd waant, maar 
waardoor men soms toch in de steek wordt BEE. We beschouwen 
de volgende situatie. 

Van een functie wenst men de nulwaarde benaderd te bepalen. 
De functiewaarden zijn bekend in een aantal punten in de buurt 
van 0. Hieruit blijkt dat de functie vrijwel lineair verloopt in deze 


1) Lezing gehouden voor de winterbijeenkomst van het Wiskundig Genootschap 
op 3 januari 1967 te Haarlem. 

De schrijver memoreert met grote. erkentelijkheid de discussies met de heer 
J. R. Zweerus tijdens de voorbereiding van deze lezing. 
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fig. 1. 


omgeving, en dat de gevraagde nulwaarde ook daarin ligt (fig. 1). 
Weliswaar zou dus lineaire interpolatie al een heel aardig resultaat 
gegeven hebben, maar men wenst nog iets beter te doen, en be- 
paalt de kwadratische functie ax? + bx + c, die op drie der gegeven 
punten dezelfde waarden aanneemt als de gegeven functie. Nu heeft 
deze kwadratische functie weliswaar twee nulpunten, maar een 
daarvan is (wegens het bijna lineaire gedrag van f) absoluut zeer 
groot, zodat men de absoluut kleinste wenst te hebben. En nu gaat 
het mis. Immers, zijn « en £ de wortels, |al < |#|, dan is 

b2 == a?. (a + B)? zm a?B?, dac = 4a?ap, dus dac < b°. Dus is 

Vb? — 4ac mb (we mogen zonder meer aannemen b > 0), zodat in 


— b + V(b? — 4ac) aanzienlijk ciĳjferverlies optreedt (dwz. men 
moet de vierkantswortel in veel meer cijfers nauwkeurig kennen 
dan men het uiteindelijk resultaat wenst). Dit cijferverlies treedt 
uiteraard niet op bij de bepaling van de absoluut grootste wortel, 
en hiermee is de remedie duidelijk: door de absoluut grootste wortel 


te delen op iel krijgt men de absoluut kleinste bevredigend. Dit 
a 


resulteert in de formule 
== DE — 2C 
NOOT b > 0 (of gn 
b 4 Vb? — 4ac b{L 4 V1 — 4acfb?} 


als men een formule wenst te hebben die, ongeacht het teken van b, 
altijd de absoluut kleinste wortel oplevert). Hierop dient men dan 
echter beslist niet vervolgens het op school geleerde voorschrift toe 
te passen dat wortels uit noemers dienen te worden verdreven, 
want dan krijgt men de oude formule weer terug. Dit geldt natuur- 
DN 
Va + vr 

4. Een ander probleem, waarvoor de oplossingswijze traditioneel 
geheel vast staat, is dat van het bepalen van het maximum van 
een functie op een segment [a,b]: men stelt de afgeleide g van f op 
en lost op g(x) — 0. 


lijk algemener voor uitdrukkingen van de gedaante 
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We zullen nu eerst laten zien hoe men gedwongen wordt stap voor 
stap dit traditionele pad te verlaten (5 t/m 7) en ons daarna con- 
centreren op een geheel ontraditionele aanpak. 

5. Zij om te beginnen g(x) —= 0 eens niet prettig oplosbaar. Dan 
kan men het numeriek proberen, bijv. met de methode van Newton 
(zie fig. 2), een inmiddels eigenlijk ook alweer traditioneel geworden 
aanpak: Uitgaande van een geschikt” punt x genereert men een 
rij {x,} als volgt: xr = Xn — 2(%)/2 (7) (dew.z. #4, is het snij- 
punt van de as met de raaklijn in het punt (%,, g(x‚)) aan de 
grafiek). Als g en 2’ voldoende netjes zijn, convergeert de rij {x,} 
snel naar het nulpunt (deze methode toegepast op g(x) = x° — a 
levert een bekend proces voor de benadering van 1/4). 





Kraj Ke 


fig. 2. 


Zodra echter de grafiek van g hobbeltjes vertoont komen de moei- 
lijkheden. 

6, Moeilijkheden komen er zeker als de afgeleide van g niet be- 
staat of niet te bepalen is, en men vraagt zich af of g(x) — O0 niet 
op te lossen is zonder g’ te gebruiken. Zo'n methode zou trouwens 
ook van belang zijn voor de nogal eens voorkomende situatie, dat 
2’ wel bestaat maar zich slechts ten koste van veel rekenwerk laat 
bepalen (de afgeleide van een functie pleegt nl. een veel ingewik- 
kelder analytische gedaante te hebben dan de functie zelf). 

Merkwaardigerwijze neemt men nu zijn toevlucht tot het in de 
analyse gewoonlijk slechts van theoretisch belang geachte proces 
van successteve halvering: zij g continu op [a, b]; zijn g(a) en g(b) 
van teken verschillend; zij m het midden van [a, b}; als g(mm) het- 
zelfde teken heeft als g(a) vervang a door m en herhaal; als g(m) 
hetzelfde teken heeft als g(b) vervang b door m en herhaal. Op deze 
wijze kan men de wortel van g(x) — 0 willekeurig goed benaderen. 

Ofschoon het successieve halveringsproces dus voor een duidelijk 
ruimere klasse van problemen bruikbaar is dan het Newtonproces, 
betekent dit allerminst dat dit laatste nu heeft afgedaan. Voor een 
groot aantal problemen werkt dit nl. aanzienlijk efficiënter. Evenzo 
treft men hiërarchieën van technieken aan om integralen te bereke- 
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nen, differentiaalvergelijkingen op te lossen etc. Steeds zullen som- 
mige van die technieken slechts vooreen vrij nauwe klasse van ge- 
vallen werken, maar dan erg efficiënt, andere minder efficiënt zijn 
maar’ algemener bruikbaar. 

7. Mogelijk echter is reeds g(= f’) niet of slechts met veel moeite 
te berekenen, zodat de vraag rijst of we het maximum van f niet 
kunnen bepalen zonder gebruik te maken van f’. Hiermee zouden 
we dan het traditionele pad geheel verlaten hebben. 

8. Onder een eentoppige functie verstaan we een niet noodzake- 
lijk continue functie die voor x S % (xQ een vast, bij de functie 
passend, getal) strikt stijgend is en voor x = x0 strikt dalend. De 
meeste in de praktijk voorkomende functies f voldoen in principe 
wel aan deze eis, althans op een geschikt interval. De moeilijkheid 
is echter dat men de functiewaarden vrijwel nooit exact in handen 
krijgt, bijv. tengevolge van afrondfouten (die men doorgaans al 
maakt zodra men de functiewaarden door een getal van een zeker 
aantal decimalen voorstelt, maar veel meer nog in de berekeningen 


fig. 3. 


die ons tot de functiewaarden voeren) en/of van waarnemingson- 
nauwkeurigheden (nl. als de functiewaarden door meting verkregen 
worden). Als men de verkregen functiewaarden grafisch uitzet komt 
er dan ook zo iets als fig. 3. Het is alsof er op de ware functie een 
soort ruisfunctie is gesuperponeerd, en de eentoppigheid is verloren 
gegaan. 

_ Bij vrijwel alle numerieke processen heeft men met dergelijke 
ruisverschijnselen te maken, en vaak veroorzaken deze aanzienlijke 
complicaties. Dit is ook het geval bij de nu te bespreken processen 
ter bepaling van het maximum. 

9. Zij f nu een op het interval [a, b] eentoppige functie waarvan 
we het maximum wensen te bepalen, waarbij we er evenwel mee 
zullen rekening houden dat de opgeleverde functiewaarden met 
ruis belast zijn. | 

Neem nu y en z willekeurig in [a, b] (zie fig. 4), y < z. Als f(z) = 
f(y) ligt het maximum op [y, 6], als f(z) < f(y) ligt het op [a, z}; 
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herhaal het proces op dit kortere interval (de situatie f(z) = f(y), 
waarbij we het interval zouden kunnen inkrimpen tot [y,‚ zl, be- 
schouwen we niet apart wegens het zeldzame voorkomen ervan). 
De vraag is nu uiteraard hoe y en z te kiezen zodat een zo klein 
mogelijk interval resulteert. Dit is kennelijk het geval voor y en z 
zo dicht mogelijk bij, en op gelijke afstanden ter weerszijden van 
het midden van [a,b], dus y =a + #(1 — e) (b — 4), 
z=ad}(l+e) (b — a), e zo klein mogelijk. Weliswaar kan een 
veel Gen interval resulteren wanneer men bijv. y en z beide 
vlak bij b legt, maar dan moet men het geluk hebben dat [y‚ 6] 
uitgekozen moet worden; heeft men pech, dan resulteert het in- 
terval [a,z], waarvan de lengte dan veel groter is dan 4(1 + €) 
(b — a). Het kleinste interval, waarvan men kan garanderen dat 
[a,b] erop te reduceren is, ongeacht de ligging van het maximum, 
heeft derhalve de lengte $(1 + €) (b — a), en wordt bereikt voor de 
gegeven waarden van y en 2; men noemt het zo ontstane proces 
e-halvering. 

Is er wel een positieve ondergrens voor e? Meestal wel, nl. t.g.v. 
de ruis. Want door de ruis kan het gebeuren dat, voor yenz dicht” 
bij elkaar, f(y) > f(z) terwijl toch het maximum rechts van z ligt; 
in dat geval zou men dus met het verkeerde interval verder gaan, 
en het maximum geheel mislopen. Men mag dus bij ruis y en z 
inderdaad niet willekeurig dicht bij elkaar nemen. Helaas kent men 
deze ondergrens voor e gewoonlijk niet, waardoor het EO moel- 
lijk hanteerbaar wordt. 

10. Nadere beschouwing van fig. 4 levert de gedachte o op dat, 
als we y en z nu toch eens niet dicht bij het midden leggen, [a, z] 
en [y, 6] wel wat langer worden, maar dat daarop dan ook al een. 
functiewaarde bekend is in een punt vrij ver van de vand, zodat men 
met het berekenen van nog slechts één functiewaarde het interval 
alweer aanmerkelijk kan verkleinen. Men krijgt dan bijv. achter- 
eenvolgens een rij intervallen [a,, b;] zoals in fig. 5, waarin we de 
optredende abscissen aanduiden met y,z,; (zodat in het geval van 
fig. 5 23 — Ya — Zj — Y, bz —.bj — z etc.), en waarin met een kruisje 
is aangeduid een absciswaarde waar f moet worden berekend. 





EE 
q | Z 0 
1 | 
; | | 
1 JV: Al b, 
a EVEN 
d, V2Z2 hb, 
EE 
449,23 b, 
fig. 5 


11. Een prettig proces, zeker voor automatische verwerking, 
zou uiteraard ontstaan als elk volgend interval (zie fig. 5) met het 
erin gelegen overgeërfde punt een verkleinde copie was van zijn 
voorganger. Dan verkeert men nl. op een schaalfactor na (die evtl. 


mmm 
a Z b 


fig. 6. 


negatief mag zijn) steeds in dezelfde situatie, zeg in de situatie van 
fig. 6, op het moment dat men een tweede punt zal kiezen. Dwz. 
dat in fig. 5 


(z — a): (b — a) = (z, — a) : (bj — 41), dus 
(2—a):(b—a)=(y —a): (2 — 4). (1) 
Als echter met [y, b] had moeten worden verdergegaan, i.p.v. met 
[a, z), hadden we analoog gekregen 
(2 —a):(b —a) == (b—z): (b — y). 
Dus (y—a):(z—a)—= (b —z): (b — y), waaruit volgt 
y—a=b—z (2) 
(welke symmetrie men ook verwachten zou). 
Dit invullen in. (1) geeft 
(2 — a): (b—a) = (b —z): (z — 4), 
de evenredigheid. voor de gulden snede. Dus 
z—a=Hb—a) (— 1/5) we 0,62 (b— a) 
y—a=th(b—a) (3 — 1/5) » 0,38 (b — a) | 
12. Bij de guldensnede-methode wordt de lengte van het inter- 
val bij elke stap met o = 4(— 1 + 1/5) ze 0,62 vermenigvuldigd, 


zodat elk vijftal stappen de lengte van het interval ruwweg tot een 
fractie reduceert. 
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De factor 0,62 kan beduidend groter zijn dan de factor $(1 + €) 
bij e-halvering, zodat men bij het guldensnedeproces mogelijk aan- 
zienlijk meer stappen zal moeten uitvoeren om een bepaalde reductie 
van het interval te krijgen dan bij e-halvering. Men zal derhalve 
zeggen dat het guldensnede-proces 1 mathematische zin langzamer 
convergeert dan e-halvering. Evenwel hoeft men bij het gulden- 
snede-proces voor elke stap (behalve voor de eerste) slechts één 
functiewaarde te berekenen tegen twee bij e-halvering. Aangezien 
het aantal te berekenen functiewaarden vrijwel steeds (nl. behalve 
wanneer f zeer triviaal te berekenen is) bepalend is voor de reken- 
tijd, kan men van het guldensnede-proces twee stappen doen in 
nagenoeg evenveel rekentijd als één bij e-halvering. Men krijgt dan 
een reductie tot een fractie a? ey 0,38, hetgeen aanzienlijk kleiner 
is dan 4(1 + e). Naar rekentijd gemeten convergeert het guldensnede 
proces dus juist sneller dan e-halvering. 

De opmerking moge triviaal zijn dat men bij praktische bereke- 
ningen, in plaats van te streven naar snelle mathematische conver- 
gentie (dus minimalisering van het aantal stappen), dient te stre- 
ven naar snelle convergentie in de tijd (dus minimalisering van het 
produkt van aantal stappen en rekentijd per stap); een feit is even- 
wel dat vele beginners deze ondekking niet aanstonds doen. 

13. In 9 constateerden we dat men i.v.m. ruis y en z niet te dicht 
bu elkaar mag nemen, op straffe van het verkeerde interval te 
kiezen. Alvorens ons af te vragen of het guldensnede-proces aan 
deze eis voldoet merken we op, dat men bij elke strategie om y, en 
z; te kiezen op den duur wel gedwongen zal zijn y, en z; te dicht bij 
elkaar te nemen, zodra nl. de lengte van het interval [a,, b,) kleiner 


hdd heten deo nch romen enke deren erker mond ontmnn omdand eerd vaansensend 
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dan e (zie 9) geworden is. Op zichzelf is dit niet verontrustend. 
Immers is tengevolge van de ruis de ligging van het maximum niet 
meer eenduidig bepaald, en is er een heel interval, waarvan alle 
punten er in gelijke mate aanspraak op kunnen maken als abscis 
van het maximum te worden aangezien; men noemt dit wel het 
onzekerheidsinterval voor bedoelde abscis. (In fig. 7 is dat het 
interval [#,q)). Het beste wat men van een strategie mag ver- 
wachten is dat hij (in de zin van de ruis) y, en z, niet te dicht bij 
elkaar neemt voordat [a,, b;] in lengte vergelijkbaar is met het on- 
zekerheidsinterval. Helaas is een zodanige strategie niet bekend, 
en ook het guldensnede-proces presteert dit niet. Zo ziet men in 
fig. 7 een situatie uitgebeeld waarin het duidelijk mis gaat. 

Tegen dit soort mislukkingen kan men zich op de in de wiskunde 
gebruikelijke wijze beschermen door het opleggen van extra voor- 
waarden aan f. Dit is hier zelfs heel eenvoudig. Theoretisch wordt 
het proces dan waterdicht. De lekkage bij de uitvoering van het 
proces houdt men echter, omdat men doorgaans niet in staat Is te 
verifiëren of de extra voorwaarden wel vervuld zijn. Een analoge 
situatie treft men aan bij vrijwel alle praktische processen. 

14. Eenmaal geraakt door de zin voor economie vraagt men zich 
vanzelf af wat in 9 nu de beste strategie is om de punten te bepalen 
waarin men f zal berekenen, d.w.z. de strategie die met een gegeven 
aantal (zeg ») te berekenen functiewaarden een zo sterk mogelijke 
reductie van het interval teweeg brengt. Deze vraag is in wezen 
dezelfde als de vraag naar de strategie die in staat is een interval 
van zo groot mogelijke lengte na berekening van # functiewaarden 
te reduceren op een interval ter lengte S.1. Dat een dergelijke zo 
groot mogelijke lengte en een bijpassende strategie bestaan, kan 
men aantonen onder de conditie dat geen twee der te kiezen punten 
minder dan een positief getal e van elkaar verschillen. 

We zullen het langs klassieke lijnen der analyse verlopende bewijs 
(compactheids-argument) hier achterwege laten. De bedoelde zo 
groot mogelijke lengte zullen we met Z, aanduiden, de bijpassende 
strategie met S 

Zij nu [a, b] (zie fig. 8) zo’n interval ter lengte /,. Laten x, en xs 
(X, > %) de eerste twee der door S, bepaalde punten zijn waar f 
moet. worden geëvalueerd. Op [a,, bj) (=[a, x,). of [xa, b]) mogen 
nu, naast de overgeërfde functiewaarde; nog # — 2.functiewaarden 
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worden bepaald. Dus x, — a Sl, b — %o Sl, (met het optre- 
den van het < teken’ moet serieus rekening worden gehouden om- 
dat de overgeërfde abscis wel eens een andere zou kunnen zijn dan 
door S,_, verlangd wordt). 

In het geval [a,, b,) = [x3, b] zijn er voor x3 twee emnerdnheder 
Ka < Kg < Xi Of 4, < %3 < b; in het eerste geval zal, op grond van 
het zojuist afgeleide, 5 — x3 Sl, In het tweede geval b — xj S 
dus steeds b — x, S be 
Wegens b —a=% —atb—x, geldt dan 

l, S a re bz (3) 
Voorts Nd (4) 

De vraag is uiteraard of in (3) het gelijkteken te realiseren is; 
in dat geval zou, gezien de afleiding van (3), moeten gelden 5 —x, = 
Xi —a=ln b— A= bg en dus ook 4, — a =l_3. We pro- 
beren het maar eens, en definiëren hiertoe een rij A, als volgt. 

A=l p= (5) 
Âr = Àra + Ârz (6) 
We nemen een interval ter lengte A, en leggen de deelpunten x, en 


Xs op afstanden A, van de uiteinden (fig. 9). Dan is 4x, — a = 
An — Ana = Âni (zie (6)) en dus x, — x3 = Ani — An-2 = Ân-3 


n—3 


Om nnn De) 
q Ân=2 Xj Xp Àn-2 b 
fig. 9. 
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Is nu [a,, b,) = [a, x‚}, dan is de lengte daarvan A, en het over- 
geërfde punt x, ligt al op afstand A3 van het uiteinde; het nog te 
kiezen deelpunt x3 op afstand A„_3 van het andere uiteinde leggen 
voert ons terug in de uitgangssituatie, zij het met alle indices 1 ver- 
laagd. We kunnen het proces dus herhalen, en eindigen na het be- 
rekenen van # functiewaarden in de situatie van fig. 10, waaruit 
inderdaad een interval ter lengte 1, waarop het maximum wordt 
aangenomen, kan worden bepaald. 
| 
4-2 bn-2 


fig. 10. 
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Ook al is het optimale proces hiermee aangegeven, prettig is het 
niet, omdat de e weer optreedt, en omdat men van tevoren » moet 
vastleggen. We vragen ons dan ook af hoeveel economischer dit 
proces zou zijn dan het zoveel prettiger guldensnede-proces. 
15. Bij de beantwoording van deze vraag maken we kennis met 
een intrigerende getallenrij, die men in uiteenlopende gebieden van 
de wiskunde tegenkomt, nl. de getallen van Fibonacci, alsmede 
met het expliciet oplossen van een recursie relatie van een algemeen 
en veelvoorkomend type. 

De getallen van Fibonacci F, worden gedefinieerd door 


Fo= Fi=lh Fe Fia + Foo k E22. 


Derhalve geldt /, < F, voor elke e. We trachten een expliciete uit- 
drukking voor de Fte vinden. Oplossingen van de recursie 


Íx = Îrar + Írm2 | (7) 


laten zich gemakkelijk aangeven: probeer maar f, — 4%, dan moet 
xk == hl J 4PT2, dus 22 = x + 1, dus, 409 = H(l + v/5) (iets der- 
gelijks kan gedaan worden voor alle homogene lineaire recursie re- 
laties met constante coëfficiënten). Nu is een lineaire combinatie 
van oplossingen van (7) weer een oplossing. Door nu c, en c, op te 
lossen uit c,‚ + ca = 1, Ci%j + C2%3 = 1 krijgt men 


rel" 


Aangezien de eerste term tussen [] snel stijgt als 2 groeit, terwijl 
de tweede dan snel daalt, geldt 


k+1 
F, en le 


rl 7 

Wegens 4, < F„ kan een proces volgens 10, waarbij » functie- 
waarden mogen worden bepaald, dus nooit garanderen dat het 
oorspronkelijke interval tot een fractie S I/F, wordt gereduceerd. 
Anderzijds geeft de guldensnede-methode in dat geval een reductie 


| 5 n—l 5 
5 . Wegens 


tot een fractie GS, =( 5 


1 /1 5\° 
GS IF Vs AM as 1,17 ziet men dat de guldensnede- 
methode welhaast optimaal genoemd mag worden, en uit 
GS al F„ * =0,62 X 1,17 ms 0,73 blijkt dat de guldensnede-met- 
hode met » + 1 punten een sterkere reductie garandeert dan welk 
proces ook volgens 9 met # punten. 
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16. Terug naar het guldensnede-proces (zie 11). De meest voor 
de hand liggende uitvoering hiervan is wellicht als olen (in het 
geval van fig. 11). | 


x= a + o(b — a) met o= H(— 1 + /5) 
x= b — o(b — a) 4 5 (8) 
X3 = % — O(4 — 4) | 

etc. 


fig. 11. 


Wie echter recht in de leer is in de theorie der gulden snede, dan 
wel zich (2) herinnert, zal zich het voortdurend vermenigvuldigen 
met d als volgt besparen: 


x= a + o(b — a) 


kp == at b— 4 (9) 
Kz = A + Xi — X3 


| 


etc. 


Dit theoretisch inzicht wordt echter niet beloond: wie het proces 
uitvoert volgens (9), ervaart dat het overgeërfde punt al spoedig 
helemaal niet meer in de buurt ligt van waar het zou behoren te 
liggen. Het komt zelfs in de verkeerde helft van het interval terecht, 
zodat elke referte aan het guldensnede-proces illusoir is. Met name 
kunnen hierdoor het overgeërfde en het met behulp daarvan be- 
paalde nieuwe punt wel willekeurig dicht bij elkaar komen. Had 
men dit kunnen voorzien? 

Aangezien o —= H(— 1 + 1/5) een irrationaal getal is, kan men 
bij het numeriek rekenen weinig anders doen dan x, = a +o,(b —4) 
te nemen met o, een goede benadering van o. Definiëren we a, door 
X3 = A + 03(%j — 4), dan volgt wegens 4 = a +b — % dat 0, = 


1 
ENE: 
01 


Analoog vinden we voor % op '[a,%] %3 = 4 + o3(% — a) met 


1 
03 = — — l (en evenzo voor 4, op [#3, %]) etc. 
Oz 1 
De recursie op; =— — l wordt fraai weergegeven in fig. 12. 
Oz | 
Het getal os wordt bij deze recursie gereproduceerd, en treedt dan 
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1 
ook op als abscis van een snijpunt van de grafieken van — — l en x. 
nj 


Men ziet dat a, zich voor stijgende k wegspoedt van o, om zich 
te begeven naar oc’ = — H(Ì + 1/5). De snelheid waarmee o, van 


1 
o wegloopt is eenvoudig te schatten door — — 1 in de omgeving 
Xx 


van o te lineariseren, d.w.z. te vervangen door de functie 


1 


— — 1 — — (# — 0), waarvan de grafiek een rechte is die in het 
0, 0 


1 1 
punt (o, — — 1) aan de grafiek van — — 1 raakt. Men ziet dan 
5 Xx 


1 
dat 0, — 0  — — (opa — 0)  — 2,6(0,1 — 0). Dus o, — 0 » 
0 


(— 2.6) (o, — 5). 


LA 
ds 


O7 g' 





EE dee er vn en DE EE OE ED EE EN en DD On EE PE ED EE DE ED 
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Hiermee kunnen we nu de omvang van het onheil voorspellen. 
Zij eens lo, — ol = 10 (heel redelijk wanneer men 1/5 m.b.v. een 
log. tafel bepaalt). Dan kan het reeds gebeuren dat 0, < #, zodat 
reeds x,, in de verkeerde helft van het interval kan liggen (dus 
zeker exit gulden snede) terwijl de lengte van het interval dan pas 
op 1 procent is gereduceerd. Wat met name verrassend is aan dit 
resultaat, is het vroege tijdstip van het deraillement. Stel nl. voor 
het gemak de lengte van [a, b] op 1. Dan verschilt x, minder dan 
10-5 van de waarde die hij zou moeten hebben. Men zou dus waar- 
schijnlijk op zijn vroegst moeilijkheden verwachten tegen de tijd 
dat de lengte van [a,, b‚] in de buurt van 10-° komt. Maar de moei- 
lijkheden komen al bij 10? Terwijl x, slechts een fout van 10 
vertoont, verschilt x,, al omstreeks 10-® van de waarde die hij zou 
moeten hebben. 

Het gesignaleerde verschijnsel is een kras staaltje van ustabilateit: 
doordat men de door de wiskunde voorgeschreven handelingen niet 
exact uitvoert (afrondfouten) kan een proces zich bij praktische 
uitvoering zeer snel zeer ver verwijderen van wat het theoretisch 
verondersteld werd te doen. Men moet op de mogelijkheid van in- 
stabiliteit steeds bedacht zijn bij het opstellen van processen, en 
het is natuurlijk de kunst stabiele processen te bedenken. Zonder 
bewijs delen we mee dat het guldensnede-proces uitgevoerd volgens 
(8) stabiel is; het gaat dan inderdaad goed totdat de lengte van 
[a;, b‚] in de buurt komt van de fout in 4, | 

17. Bij iteratieve processen rijst altijd de vraag, wanneer op te 
houden, dus naar het zgn. stopcriterium. Men kan bij het gulden- 
snede-proces bijv. ophouden wanneer het verkregen interval een 
lengte heeft die kleiner is dan een gespecificeerde waarde. Vaak 
echter is de bedoeling het antwoord zo nauwkeurig op te leveren 
als mogelijk is met het gegeven proces en de precisie waarmee de 
functiewaarden te bepalen zijn. In dat geval kan men doorgaan 
totdat een interval [a,b] verkregen wordt met daarop de punten 
Vi Zx (zie fig. 5) zodat de voor a;, vj, Zj, by bepaalde functiewaarden 
in strijd zijn met de conditie voor eentoppigheid (deze conditie 
luidt (f(z) 2 /(y,) > f(a,)] Vv (fwa) =f(ze) > f(b,)1)- Men kan ho- 
pen aldus te eindigen met een interval waarvan de lengte in grootte- 
orde gelijk is aan de onzekerheid waarmee t.g.v. de ruis de abscis 
van het maximum bepaald is; door toevallige omstandigheden kan 
het uiteindelijke interval echter wel veel korter of langer zijn (zie 13). 


Besluit: In het voorgaande zagen we hoe men ook reeds bij het 
toepassen van zeer eenvoudige wiskunde in dito processen in aan- 


158 


raking kan.komen met problemen betreffende economie, strategie, 
stabiliteit, ruis, stopcriterium. Bij vrijwel elk numeriek proces ont- 
moet men enkele van deze problemen, en het is de taak van de 
numericus tussen de klippen door te zeilen. Het zijn o.a. deze as- 
pecten, die de theoretische wiskunde niet kent, en die medebepalend 
zijn voor de identiteit van het onderhavige vakgebied; die ook een 
kenmerkende wijze van omgaan met de wiskunde dicteren. 

De vraag rijst wellicht welke rol de computer nu speelt in dit 
geheel, waar het immers duidelijk is dat de gesignaleerde proble- 
matiek gebonden is aan het proces, en optreedt zodra men dit pro- 
ces gaat uitvoeren, om het even of dit nu met de hand of met de 
computer gebeurt. 

Door de komst van de computer echter is deze problematiek pas 
actueel geworden. Een van de redenen hiervoor is ongetwijfeld. dat 
het numeriek uitwerken van problemen door deze komst zozeer is 
aangemoedigd en daardoor zoveel vaker gebeurt. Een meer fun- 
damentele reden is echter dat bijv. instabiliteitseffecten vaak eerst 
manifest worden na vele rekenstappen, veel meer dan men er met de 
hand ooit zou uitvoeren, maar veel minder dan men de computer 
wenst te laten uitvoeren. Een andere is, dat men bij het rekenen 
met de hand, al rekenende, de calamiteiten veelal ziet aankomen, 
en daar ter plaatse iets aan doet, terwijl men van de computer 
veelal slechts het eindantwoord wenst te vernemen, en de tussen- 
resultaten (die de calamiteit mogelijk zouden onthullen) onder 
water laat blijven. Het is dan ook de door de computer op grote 
schaal mogelijk geworden automatisering van het rekenen, waardoor 
de genoemde problematiek sterk aan het licht gekomen is. 
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HET ZUIVERE FUNCTIEBEGRIP 
door 
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Arnhem 


In zijn artikel „Rondom het gelijkteken’’ (Euclides 42, IX, 269 — 
277) bespreekt dr. Joh. H. Wansink de volgende vragen (,„,op- 
drachten’) 


voor welke 4: 42 + 3 = 7? 
voor welke x: 2% + 2ax J- a? —= 9? 


Wansink noteert deze vragen, prof. Freudenthal navolgend, als 
volgt: 


ol 3=7 Pl + ax Hat =|. 
Als oplossingsverzamelingen geeft hij achtereenvolgens: 


de getallenverzameling {—2; 2}, 
de „verzameling van functies” {— a — 3; — a + 3}. 


Vervolgens bespreekt Wansink: ?,A „[x? 4 2ax + a? == 9), waar- 
van hij als oplossingsverzameling noemt: de „verzameling van 
functies” {— a — 3; — a + 3}. 

Tenslotte vermeldt Wansink een moeilijkheid, nl. 

„de oplossingsverzameling van ?,A ,[x? + ax = a + 1] kan zowel 
de getallenverzameling {1}; als de verzameling {1; — a — 1} zijn"’; 
naar aanleiding daarvan introduceert hij het symbool: ? f(a) 


Alvorens hierop kritiek te geven, merk ik op: 

a. Ik deel Wansinks voorkeur voor A en Vv boven V en d als 
symbolen voor de quantificatoren „voor alle” en „er iseen’ (zie 
Wansink: Didactische Oriëntatie voor wiskundeleraren I, pag. 115); 
ik geef er de voorkeur aan, deze voorkeur metterdaad te laten blij- 
ken. 

b. Het lijkt mij van grote didactische waarde, steeds het variabi- 
liteitsgebied van een gebonden variabele te vermelden. Vaagheden 


(159) 
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kan men zich slechts veroorloven tegenover een publiek dat in het 
vak ingewijd is en dat weet wat met de vaagheden wordt bedoeld; 
een dergelijk publiek is een schoolklas maar zelden. 

Laten we nu de door Wansink genoemde vergelijkingen met de 
daarbij behorende oplossingsverzamelingen één voor één bekijken. 


biër [x? J3= 7) 
Omdat: Aserl(t? + 3= 7) ee (== ZV =—2)|, is {—2; 2} de 
oplossingsverzameling van ?erlt? + 3 = 7}. We noteren dit 
laatste zó: {— 2; 2} — {xe R|x? +3 = 7}. 


U. Prerlt® + 2ax a? — 9] 

Deze vraag heeft slechts zin, als hij als volgt wordt verstaan: 
Gegeven een reëel getal; we noemen het a; ?.p[x? + 2ax + a? = 9]. 
Nuis: Agenl(4? + 2ax Ha? = 9) ee (£= —a— Bv = —a + 3)]; 
en dus: 

{— a — 3; — a + 3} is de oplossingsverzameling van 

Prenlt* + 2ax + a? = 9); dit laatste is waar voor elke ae R,‚ en dus: 
Aaerl{#e R/x? + 2ax + a? = 9} = {—a— 3; —a + 3}. 


Ul. PrerAgerlt? + ar =atl] 


Omdat Naer (A aerl#? +ax=atll)e(4= lx? = 1), is {1} 
de oplossingsverzameling van iii. 


iv. PrerAaerlt? + 2ax J a? — 9] 
De oplossingsverzameling van iv. is leeg, omdat 
Asger (Aaerl#? + 2ax + a? = 9) ee (1 = On =OnA?= 9). 


Het ziet er naar uit, dat Wansink met de mededeling: 
{— 4 — 3; — a + 3} is de oplossingsverzameling van 
PA alt? + 2ax + a? = 9] iets anders bedoelt, dan hij feitelijk uit- 
drukt, en wel: 

Aaerl{— 4 — 3; —a +83} is de oplossingsverzameling van 
Perl? + 2ax + a? = 9]; en dat hij met de bewering: 
{1; — a — 1} is de oplossingsverzameling van ?,A „[x2 + ax =a + 1, 
feitelijk bedoelt: 

A aerl{l; — a — 1} is de oplossingsverzameling van 

boerlt +44 =a + 1]. 


Het is maar een kwestie van de plaats van de quantificator! 


Vv. De „verzamelingen van functies’’ en het symbool ?…_ va) 


Op pagina 273 schrijft Wansink: 
Terwijl de oplossingsverzameling van de vergelijking ?,(x2 + 3 == 7) 
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een getallenverzameling was, beschouwen we een verzameling 
van functies van a als oplossingsverzameling van de vergelijking 
(2 + ax Hat=9). 

Substitutie van — a — 3 en van — a + 3 doet het oordeelsschema 
x° 4 2ax + a? = 9 overgaan in een identiteit van a. 

Wansink suggereert dat — a — 3 en — a + 3 functies van a zijn, 
en dat er twee „functies x van aq’ zijn die 4x? + 2ax + a? = 9 
doen overgaan in een identiteit. 

Laten we de zaken nauwkeurig bekijken. 

We noemen de klasse van de functies van R naar R: F. 

Wansink schijnt met zijn vraag op pagina 274: 

Peep) [A° + Zar a* = 9] te bedoelen: 


toerAaerlP?(4) a Za. pla) Er a? EE 9), (1) 
Nu is: 
A ger [A aer(P°(4) + 2a. pla) Ja? = 9] > 
A ger (P(4)= —4a—8V pla) ie + 3]; 
en dus- | 
pe FIA serlp?(a) + 2a. pla) + a? = 9} = 
we FIA gerlpla) = — a — 3V pla) = — a + 3). 


De oplossingsverzameling van (1) bevat dus oneindig veel elementen, 
hieronder zijn natuurlijk de twee functies die Wansink waar- 
schijnlijk bedoelt, nl. x — — x — 3 en Xx — — X + 3 maar ook bij- 
voorbeeld de functie g : R — R die gedefinieerd wordt door: 


pi (x) = — x — 3, als x rationaal is 
g(x) = — x + 3, als x irrationaal is. 


Tot besluit: | 
a. Ik twijfel er niet aan, of men kan leerlingen in korte tijd leren 
oplossen: 


xe taxtl=0et=— dat Va— dvi Ja $Var—d. 


_ Men geeft hun daarmee geen idee van wiskundige precisie. 

Ik denk, dat het Wansinks bedoeling was, met zijn symbool 
 e=f(a) En met zijn opmerking ‚dat ook moet worden opgegeven voor 
welke deelverzameling van R de te bepalen functies zijn gedefini- 
eerd (in dit geval de deelverzameling {ae R|a S —2vaz?2)’, 
die precisie enigszins te herstellen. 

Het is jammer, dat bij zijn voorstel het zuivere functiebegrip in de 
mist verdwijnt. Ik kan zijn voorkeur voor het symbool hen f(a) 
dan ook niet delen. 
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b. De zondebokken zijn hier zonder twijfel opdrachten zoals: 
„Gegeven x? + ax + 1 = 0; druk x uit in a” (ofwel: los x op uit: 
x° + ax + 1 = 0) en 
„Gegeven x? + ax + 1 = 0; druk a uit in x” (ofwel: in a op uit: 
x2 + ax J- 1 = 0). 
Hiermee kan -dunkt me- niets anders bedoeld zijn dan: 
„Zij ae R; voor welke zeR: x? + ax + 1 = 0?" resp. 
„Zij xe R; voor welke aeR: 4? + ax 4 1= 0?" 
De antwoorden luiden: 
Als — 2 <a < 2dan: {xe R|x? + ax + 1 = 0} =g, 
en als a< — 2va 2 dan: weR|x? + ax +1 = 0} 
{hat HVA 4 — Ja — Har 4, 
respectievelijk: 
Als x = 0, dan: {ae R|x? + ax + 1 =0}=@, en « 





— Xx —1 
voor alle x waarvoor x # 0: {ae R|x? + ax + 1 = 0} = 5 
% 


Op deze wijze zijn de problemen helder gesteld en worden eerlijke 
oplossingen van die problemen gegeven. Hierbij speelt het functie- 
begrip vooralsnog geen rol. Wie dat wel een rol toebedeelt, wake 
er voor, niet in het moeras te geraken dat onder v. is gesignaleerd; 
ik vrees, dat middelbare scholieren in behendigheid en kracht te 
kort schieten, om zich uit dât moeras te redden. | 

De opdrachten „druk x uit in a” resp. „druk a uit in #” zijn te 
wazig om in-het-begin-toelaatbaar te kunnen worden geacht: we 
beschikken immers over een heldere formulering. 

Om natuur- en scheikundecollega’s terwille te zijn, kunnen wis- 
kundeleraren deze vage uitdrukkingen natuurlijk introduceren, 
wanneer eenmaal dergelijke problemen in hun heldere vorm door 
de leerlingen begrepen worden. Het is goed, dat leerlingen ervaren 
dat „„ook wiskundigen slordig zijn, maar dat wiskundigen tenminste 
wéten, dat ze het zijn’’ (naar een uitspraak van prof. van der Blij). 


Naschrift. 


(1) Ik ben de heer Maassen zeer dankbaar voor zijn indringende, 
kritische analyse van het onderhavige probleem. Zijn bezwaren zijn 
van logisch standpunt beschouwd stellig steekhoudend. Ik hoop, 
dat het inderdaad mogelijk zal blijken. onze leerlingen te brengen 
tot het oplossen van vergelijkingen met parameters volgens metho- 
den die leiden tot de-door Maassen voorgestelde precisie van uit- 
drukking. Zijn doel zal echter naar mijn mening slechts bereikt 
kunnen worden, als alle vergelijkingen met parameters aanmerke- 
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lijk later aan de orde worden gesteld dan thans gewoonlijk het geval 
18, 
(2) Ik aanvaard Maassens bezwaren tegen de plaats van de kwan- 


toren geheel. Zie voor mijn opvattingen de opmerking op p. 265 
van deel II van mijn Didactische Oriëntatie. 


(3) Het verschil in opvatting tussen Maassen en mij voor wat 
betreft het symbool x = f(#) onder de interrogatieve kwantor hangt 
samen met de omstandigheid, dat mij hierbij nog voor ogen heeft 
gestaan het functiebegrip geformuleerd in het rapport van de No- 
menclatuurcommissie 1959 op p. 55 van Euclides 35. Hier lezen we 
namelijk nog: „Men noemt 2x + 5 een functie van x'’. In het wis- 
kundig onderwijs van vandaag en morgen komt voor dit primitieve 
functiebegrip stellig het door Maassen gehanteerde in de plaats. 


(4) Op grond van de geopperde bezwaren zou ik thans het sym- 
bool x = f(b) onder de interrogatieve kwantor liever door x(#) 
zien vervangen, ook al begrijp ik, dat ik daarmee aan Maassens 
wezenlijke bezwaren niet tegemoet kom. Het symbool krijgt dan de 
‘ betekenis van „druk x in p uit”. 
Ik vind de opdracht x in p uit te drukken (en voor de uitvoering 
van deze opdracht staat men ook bij Maassens methode) van een 
voor ons beginonderwijs acceptabele vaagheid. Ideale precisie kome 
aan het slot, men starte er niet mee. Het eigenlijke functiebegrip 
blijft hierbij dan. nog enigszins op de achtergrond. Er wordt echter 
aan een latere precisering in Maassens geest niets in de weg gelegd. 
Ook in de techniek van het oplossen van vergelijkingen met para- 
meters in de geest van Maassen is het vinden van „uitdrukkingen”’ 
als — Ja + Ay/(a? — 4) en — ha — by/(a? — 4) een wezenlijk 
element. Dat daarbij het zuivere functiebegrip tijdelijk nog in de 
mist blijft hangen, vind ik niet catastrofaal, mits de docent maar 
zorg draagt, dat die mist voor de leerling te zijner tijd optrekt. 


Joh. H. Wansink 


VERSCHEIDENHEDEN 
door 
Prof. dr. O. BOTTEMA 
Delft 


LXX. Elementary, dear Watson. 


De thans zo onaanzienlijke geometria, assepoester in het huidige 
mathematisch bestel, krijgt toch nog op haar tijd fanmarl uit de 
smalle gemeente van haar bewonderaars. De eenvoudige plani- 
metrische opgave, die Frans van Schooten meer dan drie eeuwen 
her aan Christiaan Huygens voorlegde en die in Verscheiden- 
heden LXVII (Euchides, 42, 1966/67, 204— 208) werd besproken, 
bleek aanleiding te geven tot menig schriftelijk en mondeling weer- 
woord. Met toestemming van de redactie en van de hieronder met 
name genoemde correspondenten worden in deze bijdrage de op- 
merkingen samenvattend verzameld. 

Het vraagstuk luidde als volgt. Gegeven is de driehoek A BC met 
op AB en AC respectievelijk de punten Den E. Gevraagd wordt op 
BC een punt F te construeren, zo dat de hoeken BDF en FEC 
gelijk zijn. Ter aangehaalde plaatse werd een drietal oplossingen 
aangegeven. In de eerste twee werd de meetkundige plaats bepaald 
van de punten S, waarvoor de hoeken BDS en SEC gelijk zijn en 
wel respectievelijk projectief-meetkundig en analytisch. Het re- 
sultaat was een orthogonale hyperbool, waarvan dan het snijpunt 
met BC het gevraagde punt F is. Een derde methode bracht de 
vraag terug tot het oplossen van een goniometrische vergelijking 
van een bekend type. 

Ten aanzien van de mogelijkheid van een elementaire oplossing was 
het in de gedachtengang van de schrijver niet verder gekomen dan 
tot een vaag vermoeden. 

Wij geven nu maar dadelijk de oplossing van Dr. J. T. Groen- 
man, die naar de mening van alle betrokkenen de prijs zou mogen 
krijgen voor de beste inzending. Zijn constructie luidt, na het aan- 
brengen van een kleine retouche door Prof. dr. S. C. van Veen, als 
volgt (fig. 1). 

Spiegel de driehoek in de basis BC; zij E’ het beeldpunt van ZE en 
G het snijpunt van CE’ met het verlengde van AB. Daar / FEC 
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= ZFE'C zijn /GE'F en /GDF elkaars supplement, waaruit 
volgt dat GE’ FD een koordenvierhoek is. Het gevraagde punt F 
wordt dus gevonden als snijpunt van BC met de cirkel door de drie 
bekende punten D, Gen E'. 





fig.1. 


B ligt tussen G en D en dus binnen de cirkel, C ligt buiten GE' 
en dus buiten de cirkel: tussen B en C ligt dus steeds één enkel 
punt van de cirkel. Het tweede snijpunt van de lijn BC met de cirkel 
levert de figuur die ook reeds in het oorspronkelijke opstel genoemd 
werd. | 

De oplossing laat aan eenvoud niet meer te wensen over. Dat 
rechtvaardigt als begeleiding de klassieke opmerking van het genie 
uit Baker Street tot zijn goede vriend, die een, neiging had te veel 
achter de dingen te zoeken. 

Ook Prof. dr. G. R. Veldkamp zond ons een elementaire oplos- 
sing van het vraagstuk. Daarin werd eveneens een spiegeling toe- 
gepast en wel ten opzichte van de bissectrice van 4. In eenvoud lag 
zij bij de bovenstaande constructie achter en daar het goede gaarne 
wijkt voor het betere, wordt zij hier niet gereproduceerd. Het- 
zelfde geldt voor een tweetal van M. J. Ritter ontvangen oplossin- 
gen, die het gevraagde punt door middel van een algebraïsche ana- 
lyse bepalen. 

De enige correspondent die het vraagstuk bleek te kennen was 
Drs. R. Kooistra, die zich reeds enige jaren geleden er mee heeft 
beziggehouden daarbij zowel de goniometrische oplossing als de 
analytische (met een wat ander coördinatenstelsel) vond. 
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Een uitvoerig schrijven over het onderwerp ontvingen wij van 
Prof. dr. S. C. van Veen, wiens jarenlange omgang met Huygens’ 
werk welbekend is. Hij deelt allereerst mede dat ook hem uit de 
in de Oeuvres gepubliceerde briefwisseling niet is gebleken, dat ooit 
op het gestelde probleem is teruggekomen. Zich verder met aandacht 
verdiepend in de vraag hoe Huygens het eventueel had kunnen 
oplossen, komt van Veen op grond van de denkwijze van tijd en 
milieu en gegeven Christiaan'’s kennis van zaken omstreeks 1648, 
tot de uitspraak dat de meetkundige oplossing met behulp van de 
hyperbool het meest waarschijnlijk is. De negentiende-eeuwse 
voortbrenging daarvan door projectieve stralenbundels zal hij niet 
hebben benut, maar Huygens kende zijn Apollonius en moet 
op grond van diens methodieken en theorema's de meetkundige 
plaats hebben kunnen afleiden en ook wel de snijpunten daarvan 
met de basis van de driehoek. 

In de correspondentie van Huygens komen nog vele andere 
meetkundige opgaven, van uiteenlopende aantrekkelijkheid, aan 
de orde. Het blijkt niet nodig ze uit de Oeuvvres bijeen te zoeken. Het 
is ons namelijk achteraf gebleken dat zulks reeds is gebeurd door 
P. van Geer, die onder de titel Hugenia geometrica een, twaalftal 
artikelen publiceerde over het meetkundige werk van Huygens in 
opvolgende jaargangen van het Nieuw Archief voor Wiskunde van 
1907 tot 1913. In het eerste artikel (N.A. v. W., 2e reeks, deel VII, 
1907, 215-—226) treft men uit de correspondentie van 1646 tot 
1656 in chronologische volgorde 22 dezer opgaven aan; het door ons 
besprokene vindt men op bladzijde 220 onder no. 7. Als de cor- 
respondentie ook de oplossing bevat wordt deze door van Geer 
bijgevoegd. Bij no. 7 (en een aantal andere) komt geen bijschrift 
voor, zodat het wel zeer waarschijnlijk wordt dat Huygens en zijn 
wetenschappelijke vrienden niet op de opgave zijn teruggekomen. 

In zijn algemene beschouwingen herinnert van Geer er aan dat 
noch de coördinatenleer van Descartes, noch de ontwikkeling van 
de trigonometrie door Snellius op de denkwijze van Huygens 
van grote invloed zijn geweest en dat deze verknocht bleef aan de 
zuiver meetkundige methode der ouden, een opmerking die het 
door van Veen geuite vermoeden versterkt. 


DE PARAMETERVOORSTELLING VAN EEN PUNT OP DE 
PARABOOL EN OP DE LIJN 


door 


R. KooIsTRA 
Ede 


1. In de wiskunde l.o.-rubriek van het Nieuw Tijdschrift voor 
Wiskunde en ook in de klas gebruiken we steeds de handige 
parametervoorstelling P(2pA?; 2pA) voor een punt P, gelegen op de 
parabool y? —= 2p4. Is de parameter p van de parabool bekend, 
dan is vaak P nog eenvoudiger voor te stellen b.v. 

ye =—= 8x — P(2Â2; 4); y: = —8x > P(—2Â?; 44), 
ya == 4x — P(A; 24); yi == — 4x — P(— A2; 2Â). 

In twee van de ons ter beschikking staande v.h.m.o.-leerboeken 
vonden we deze voorstelling van P vermeld !). 

Het moet de klas daarbij goed duidelijk zijn, dat het teken van 
x= 2pÂ* voor zowel p > 0 als p < 0 zonder meer in orde is 
(b>Oeex>0; Pb <Oeex <0) en dat we, omdat A toch elke 
waarde aanneemt, van de oplossingen y = + 2pÂ van y? — 4p2}? 
die met het positief-teken kunnen kiezen. | 

De voordelen van deze voorstelling zijn: 

a. Eris één variabele A en niet twee, die bas als we P(%, 9) 
zouden stellen. 

b. De betrekking #? = 2p4 kan vervallen en } kan dus ook niet 
worden an 


2. a. We nemen als voorbeeld het volgende vraagstuk: Gegeven 
vs de parabool y? = Ax. Een punt P doorloopt de parabool. De raak- 
lijn in P smjdt de x-as in A en de y-as in B. Bepaal de verzameling 
van de middens M van AB. 

le Opl.: Stel P(A? 24). De vergelijking van de raaklijn in P 
luidt Ay = x + A2, zodat de snijpunten A en B opv. zijn: 

(— Â2;0) en (O; A) (A #0). Is nu M (£#;#), dan geldt: 


X En 
=H Tjerk 


2) D.J. E. Schrek, Beknopte An. Meetkunde, blz. 87; P. Wijdenes, An. Meet- 
kunde blz. 44. 
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zodat de gevraagde vergelijking y? = — 4% is. 
2e Opl.: (zonder parameter-voorstelling). Stel P(%;#), dan is 
j2 —= 4x. De raaklijnvergelijking is nu: #y — 2x + 2%, waaruit de 


snijpunten A(— #; 0) en Blo). Is M(x; y), dan geldt: 


=d | z2 z2 
X syteem ta (2 #0, #0). 
ji 

Je 

Onderzoek van P(O; 0) en A = Oleert, dat ook de oorsprong tot de 
verzameling y? = — Jx behoort. 

Men vergelijke beide oplossingen en oordele zelf. Overigens is ook 
de 2e opl. nog vrij wat eenvoudiger dan die in het leerboek !), 
waaraan we dit vraagstuk ontlenen. 

b. Een wel heel mooie toepassing biedt ons het 3e vraagstuk van 
het eindexamen 1963: 

Bewijs dat elke raaklijn van de parabool 7 == 8x behoort tot het 
stelsel lijnen x — py + 2p? = 0. (p varabel). 

Opl.: De raaklijn in P(24?; 44) van de parabool is de lijn met 
vergelijking 


Ady — Ax + 8/2 
of 
x— Ây + 22 = 0 q.e.d. 


c. We ontlenen aan P. Wijdenes, Beknopte anal. meetkunde 
nog het volgende voorbeeld (blz. 64, vb 3): 

Voor de loodrecht op elkaar staande koorden OA en OB van de 
parabool yr —= Wx geldt: OA =8 XxX OB. Bereken de richtings- 
coëfficiënt van AB. i 

le Opl.: (zonder parametervoorstelling). Zij m de richtingsco- 
efficiënt van OA, dan is die van OB gelijk aan — 1/m. Om de ge- 
dachten te bepalen zullen we m >> 0 onderstellen; spiegeling in de 
x-as levert dan tenslotte een tweede oplossing. 

De abscis van A is de van nul verschillende wortel van (mx)? = 2px 
dus gelijk aan 2p/m?, die van B is 2pm?,- Voor de ordinaten he 
geldt: 

ap? 


2E een 2 
nn opv. YB = 4p*m?, 


1) Drs. P. E. Lepoeter, Gids voor de an. meetkunde; blz. 117 e.v. 
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Nu is: 
8 Ap. 4D? 

OAt == Et Ae 

en 


OB? = 2 J yà — 4p2 mt J Ap2m?. 
Uit OA? = 64 OB? volgt nu: 


a (1 4 m2) — 64 m?(l + m?) 


of 
me = tea > M= Ig. 
De punten A en B zijn dus opv. (85; 4p) en ($p; — p), dus is de 
gevraagde richtingscoëfficiënt: 
id 
Spip Tb 
De gevraagde coëfficiënten zijn dus + 2/. 
2e Opl.: (met parametervoorstelling). We stellen de punten A en 
B opv. (2pÂ°; 2p4) en (2pt?; 2pt). De richtingscoëfficiënt van AB is 
dan: EE 


2 
5 





De WA) 1 Ar 
Tp ari 
_ Uit m, m, = — 1 voor de richtingscoëfficiënten m, — 1/A van OA 
en m, = If van OB volgt: IM = — 1, £= — 1/A. Nu geldt: 
OA? =640B? dus  4p2(A2 4 A4) = 64. 4PU(L2 + 1) 
of: | 
(1 A2) = 64. (lS) of AB(L 4 42) = 64. 1/A2(1 + 1/42) 
waaruit: 
| 1 2 
nd 3 


3. Niet altijd verdient de parametervoorstelling de voorkeur 
boven de gewone voorstelling P(£; 9) met 7? —= 2p%. Het enigszins 
beruchte examenjaar 1965 legde de examenkandidaten de volgen- 
de opgave voor: : 

Gegeven de parabool y*= 4x en de cirkel (x — 2}? + y= 8. 
Stel de vergelijking op van de raaklijnen aan de parabool, die tot 
(4:0) de afstand 31, hebben. 
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Opl.: De parametervoorstelling P(A? 24) voert tot de vergelijking 
4 JJ? 
VI+ A? 
vergelijking overigens eenvoudig op te lossen is: A= + 1/3 of 


A= dt dV5. 
Stellen we echter P(X; #), dan is de raaklijn-vergelijking 


== 7|,, die menig leerling onrustig maakt, hoewel de 








Vy — 2x — 2X= 0 


zodat we de vergelijking 




















BEEST or LERS ot [AES] 
Vr 4 Ô VAR + 4 ì VE +1 î 





moeten oplossen, die er iets eenvoudiger uitziet: 2 = 3 of # = 5/,. 


4. Bij de afleiding van de vergelijking van de raaklijn aan de 
parabool 4? == 2px neemt menig schrijver zijn toevlucht tot de 
differentiaal-rekening, om ingewikkeld rekenwerk te ontgaan. Het 
blijkt, dat we, gebruik makende van de parametervoorstelling, 
zeer wel de differentiaalrekening kunnen missen. Onder 2 vonden 





1 
we voor de r.c. van de sniĳlijn AB: m = EE zodat zijn ver- 


gelijking luidt: 
1 
At 


Deze vergelijking gaat nu over in die van de raaklijn in A als £ — A, 
dus 





y— UPÀ = (x — 2pA2). 


1 


Stellen we nu weer A (#9) en substitueren we Â = z/y, dan ver- 
krijgen we: 


y= Je (x — £) 
of 
2pry — 22 — PVA — DYL 
Fy = PI + DE, 
zodat de vergelijking van de raaklijn in A (#;) luidt: 
jy =— bx + PE. 
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5. Ook de afleiding van verschillende eigenschappen van de 
parabool verloopt vlot, als we van de parameter-voorstelling ge- 
bruik maken: 

a) Is AB een koorde van de parabool y? = 2px, dan is de rc. 


Ì 
van de drager van AB (zie 2.): m — per: Het midden M van AB 


is M (p(A2 +42); P(A +4), waarvan de ordinaat constant is, 
als dat met » het geval is, dus: 

De verzameling van de middens van evenwijdige parabool-koorden 
vs een halve lijn evenwijdig aan de as van de parabool. 

b. De vergelijking van de raaklijn in A(2p4?; 2pA) luidt 


Uy == x J 2PÂ2, 


Het snijpunt van de raaklijn en de y-as is (0; pA), het heeft dus de 
halve ordinaat van A tot ordinaat, bijgevolg is: O het midden van PQ 
waarbij Q (2pA%; 0) de projectie is van A op de x-as en P(— 2p4?; 0) 
het snijpunt is van de raaklijn en de x-as. 

De normaal in A heeft tot vergelijking 


y— 2pÀ= — 2Alx — 24°) 


en snijdt dus de x-as in R(p + 2pA2; 0), dus de subnormaal QR —= p 
en F(4p; 0) is het midden van PR. 

c) Voor een koorde AB, die door F gaat, is de volgende betrekking 
af te leiden: 

De drager van AB heeft de vergelijking 





pW je — 2) 


AB gaat door F, als voldaan is aan 


— WPA HE) = Ib — WA 


of aan 
M= — la 
dus | 
ZEE 
Ya YB 


Onderstellen we p > O en y‚ > 0, dan is 4, > 0, zg > Oenyg < 0, 
zodat voor een koorde AB door F de betrekking geldt: 


tg AOF . tg BOF — 4. 


onafhankelijk van #. 
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6. a. Zonder nu de parametervoorstelling van een punt P (r cos p; 
ysinp) op de cirkel 42 Jy? —=7° waarbij p= / (OP; x-as), 
te willen propageren, bewijst hij toch wel eens goede diensten als 
bv. in het vraagstuk 3 van het examen h.b.s. 1962 (analytische 
meetkunde). De gewone eliminatie van x, en y, is nogal lastig: 

Gegeven zijn de cirkel x?% + yt =r? en het punt A (a;0); a > 0, 
a Hr. P(x‚, y‚) is een punt van de cirkel. Bepaal de verzameling van 
het snijpunt S van OP en de as van AP, als. P de corkel doorloopt. 

Opl.: Stel x, = 7 cos p en y, = 7 Sin p. 

De r.c. en het midden van AP zijn opv. 


EE ren 
1 COS p —- q Dn {#(r cos p SL a); r sin p}. 
| jp ke (1) 
S is het snijpunt van | (rcosp — a)x + (rsinp)y = F(r2—a?) (2) 
Uit (Il) volgt: co 4 y zie 
Ì i COSP= == eN SNP di 
Vr? + y2 Vr? + y? 


gezet in (2) geeft na vermenigvuldiging van beide leden van (2) 
met Vx? + y? en na kwadratering daarvan: 


A Hy!) = (pt — Hat +} an)?, 
een ellips of hyperbool voorstellende al naar gelang a? < r? of 
at > r2 is. 
We merken nog op, dat uit de bijbehorende figuur van deze opgave 
eenvoudig valt af te lezen: 


SO + SA =r als A binnen de cirkel ligt en 
SO — SA =r als A buiten de cirkel ligt. 


b. Ook voor een punt op een rechte lijn heeft het vaak voordeel 
de parametervoorstelling te hanteren bv. | 
P op de lijn 4 — y + 2= 0 — P(A;A + 2) (eindex. v.h.m.o. 
1966) 
Pop delijn x — 2y 4 4 —=0 — P(24 —4;Â) (eindex. v.h.m.o. 
| 1963) 
P op de richtlijn van y? == 2px — P(— Jp; À) 
(eindex. v.h.m.o. 1962) 
P op de lijn 2x + 9y —= 18 — P(OA; 2 — 24) | 
(Euclides, mrt '67, blz. 188) 
c. Tenslotte nog de parametervoorstelling P(acosp; b sin p) 
voor een punt van een ellips, die te voorschijn komt in het volgende 
vraagstuk: 
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Op een lijn l liggen de vaste punten A, P en B in deze volgorde. 
Wat is de verzameling van P als Aen B langs de benen van een 
rechte hoek glijden? ( Proclos, 410—485). 

Kiezen we het been, waarop A ligt, als x-as, het andere als y-as 

en stellen we PA = 5, PB = a en ZOAB = p, dan volgt de ver- 
zameling, zijnde een ellips, eenvoudig door eliminatie van p tussen 
Xj = A COS Pp en y, = b sin p. 


GEORG FERDINAND LUDWIG PHILIPP CANTOR 


Cantor, de schepper der verzamelingsleer, is 50 jaar geleden, op 
8 januari 1918, te Halle overleden. Op 3 maart 1845 is hij te Sint- 
Petersburg geboren, waar zijn uit Kopenhagen afkomstige vader 
zich gevestigd had. In 1856 verhuisde het gezin naar Frankfurt am 
Main. Nadat Georg op wens van zijn vader eerst twee jaar techniek 
gestudeerd had, kreeg hij toch in 1862 diens toestemming om wis- 
kunde te studeren. Eerst was hij een jaar te Zürich, daarna in Berlijn, 
met nog een korte onderbreking te Göttingen. Na zijn promotie in 
1867 werd hij achtereenvolgens privaatdocent (1869), buitenge- 
woon- (1872) en gewoon hoogleraar (1879) te Halle. 

Van 1872 dateert een kort artikel Über die Ausdehmung eines 
Satzes aus der Theorie der trigonometrischen Reihen (Math. Annalen). 
Hierin vinden we de invoering van het reële getal als limiet van een 
convergente rij van rationale getallen: de „fundamentaalrij” van 
Cantor, alsook al de kiem van zijn latere verzamelingsleer. Zijn 
grote werk daarover is in zes afleveringen in de Mathematische An- 
nalen verschenen (1879—1884): Über unendliche lineare Punktman- 
mgfaltigkeiten. Al eerder had hij in een drietal korte artikelen de 
begrippen „„machtigheid’ van een verzameling en „aftelbare”’ ver- 
zameling ingevoerd en bewezen, dat de verzameling der reële getal- 
len niet aftelbaar is en dat de machtigheid van het continuum onaf- 
hankelijk is van zijn dimensie. In het genoemde grotere werk wordt 
dit alles nog eens herhaald; verder vinden we erin de thans zo be- 
kende begrippen „doorsnede”’ en „vereniging’ van, verzamelingen, 
dan nog de definitie van geordende verzameling en ook een construc- 
tieprincipe voor oneindige verzamelingen met opeenvolgende mach- 
tigheden, met als eerste de aftelbaaär-oneindige verzameling, dan een 
met de tweede machtigheid, de derde enz. Aan het slot spreekt 
Cantor het vermoeden uit, dat het continuum de tweede machtig- 
heid heeft. Van dit laatste stelt hij een bewijs in het vooruitzicht; 
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hij heeft het echter niet kunnen geven en tot op heden is dit zoge- 
naamde „continuumprobleem’’ onbewezen. 

Cantor heeft op zijn ideeën, het invoeren van het „actueel-onein- 
dige” in de wiskunde, veel kritiek gehad, met name van zijn vroe- 
gere leermeester Kronecker, die te Berlijn een invloedrijke positie 
bezat. Dit en het feit, dat hij er niet in slaagde zijn voor de machtig- 
heid van het continuum uitgesproken vermoeden te bewijzen maak- 
ten hem neerslachtig, waardoor hij enige tijd een afkeer van zijn 
wiskundig werk kreeg. Anderzijds had hij ook wel aanmoediging en 
hij vond veel steun in de van 1872 daterende vriendschap met 
Dedekind. 

Later was hij een dergenen, die ijverden voor het tot stand komen 
der „Deutschen Mathematikervereinigung’’, waarvan hij van 1890 
tot 1893 de eerste voorzitter geweest is. 

In zijn laatste werk Beiträge zur Begründigung der transfiniten 
Mengenlehre (Math. Annalen, 1895, 1897) geeft hij nog eens een over- 
zicht van zijn verzamelingsleer, Voor machtigheid vinden we nu de 
term „kardinaalgetal’’; het eerste transfiniete kardinaalgetal, dat 
der aftelbaar-oneindige verzameling, heet „Alef-nul'’, enz. 

Cantor’s ideeën hebben tenslotte ruime erkenning gekregen, te 
te beginnen op het eerste internationale congres van wiskundigen 
te Zürich in 1897. Ook ontving hij verschillende wetenschappelijke 
onderscheidingen. De verzamelingsleer werd een aparte tak der wis- 
kunde waarvan in 1906 (W. H. Young) en in 1914 (F. Hausdorff) 
de eerste handboeken verschenen. 

A. J. E. M. Smeur 


RECREATIE 


Nieuwe opgaven met oplossingen en correspondentie over deze rubriek gelieve 
men te zenden aan Dr. P. G. J. Vredenduin, Kneppelhoutweg 12, Oosterbeek. 


189. 12-34 J- 6789 = 67-89 + 1234. 
(«) (v) (2) (#) 
Vul in zo, dat een analoog resultaat ontstaat: 
33:82 + .... = ..".. + 3382. 


Stel een methode op om bij gegeven u en v een oplossing voor # en y te vinden. 
(B. Kootstra) 


190. Gegeven is een zeshoek met daarbinnen vier punten. Verbind de hoekpunten 
van de zeshoek en de erbinnen gelegen vier punten zo door middel van lijnstukken, 
dat de zeshoek verdeeld wordt in driehoeken, waarvan de hoekpunten uitsluitend 
in deze tien punten liggen. Met hoeveel verbindingslijnstukken is dit minimaal mo- 
gelijk ? | 
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191. Een stad bestaat uit congruente vierkante huizenblokken met straten ertus- 
sen. De breedte van de straten wordt verwaarloosd. Drie winkels bevinden zich op 
drie straathoeken, zoals in de linker figuur is aangegeven. Op een straathoek moet 
men een magazijn inrichten, van waaruit de winkels bevoorraad worden. Een eco- 
nomische eis is deze straathoek zo te kiezen, dat de som van de afstanden langs de 
straten gemeten van het magazijn tot de drie winkels, minimaal is. Waar kiest men 
het magazijn? 

Zelfde vraag voor de rechter figuur, waarin vier winkels bevoorraad moeten wor- 
den. N 

Ten slotte wordt gevraagd het probleem algemeen op te lossen voor een willekeu- 
rig aantal winkels. 





(Ontleend aan een voordracht van G. Hall in Knokke, 1967) 


' OPLOSSINGEN 


187. Een tegelvloer is 8 tegels breed en # tegels lang. Een diagonaal heeft met het 
zevendedeel van de tegels een lijnstuk gemeen. Bereken z. 

Het is duidelijk, dat # = 7. De diagonaal zal 7 verticale en 7n — 1 horizontale 
lijnen snijden. Zijn de snijpunten alle verschillend, dan loopt hij dus over 7 + 7u — 1 
+ 1 tegels. Als de diagonaal door p hoekpunten gaat, wordt dit aantal gereduceerd 
tot 7 + 7 — Pp. 


Er zijn nu de volgende mogelijkheden: 


n is deelbaar door 2, maar niet door 4, en p = 1, 
n is deelbaar door 4, maar niet door 8, en p = 3, 
n is deelbaar door 8 en p = 7, 
n is niet deelbaar door 2 en p = 0. 


Deze mogelijkheden leveren resp. 


74 In — l=8nen dus” = 6, 
7 + In — 3 == 8n en dus n —= 4, 
7 + In — 7 == 8nen dus” = 0, 
7 + In —= 8n en dus n =—= 7. 


We vinden dus drie mogelijkheden, t.w. n= 6, n = 4, n= 1, en dus x = 42, 
== 28, 4 —= 49, 
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188. A neemt een getal van vier verschillende cijfers in gedachten; het mag met 
een 0 beginnen. 


‚B raadt: 0123; A zegt: 1 goed 0 raak. 
B raadt: 4567; A zegt: 1 goed 1 raak. 


Dit betekent: van de cijfers O, 1, 2, 3 komt er één in het getal van A voor en dit 
staat in 0123 niet op dezelfde plaats als in het getal van A. Van de cijfers 4, 5, 6, 7 
komt er één voor in het getal van A en dit cijfer staat op de goede plaats. 

B mag nog drie keer raden en moet dan het getal van A weten. Hoe gaat hij te 
werk? 

1. B raadt: 4523. Al naar het antwoord, dat hij krijgt, weet hij nu: 


één van de cijfers 4 en 5is raak - of - één van de cijfers 6 en 7 is raak, 
één van de cijfers 2 en 3 is goed (maar niet raak) - of - één van de cijfers 0 en 1 is 
goed. 


Bovendien weet B nu, dat de cijfers 8 en 9 in het getal van A voorkomen. 
Er zijn twee principieel verschillende gevallen: 


a. 4 of 5 raak en 0 of 1 goed, 
b. 4 of 5 raak en 2 of 3 goed. 


We beschouwen eerst geval a. 

2. B raadt: 4019. Is het antwoord 3 goed, dan was 4 raak, en is het antwoord 2 
goed, dan was 5 raak. 

Onderstel het antwoord is 3 goed 1 raak. Het getal van A kan dan geweest zijn: 


4908 4980 4890 4981 4891. 


B raadt 4908. De antwoorden, die hij in de vijf mogelijke gevallen krijgt, zijn alle 
verschillend. Hij kan dus het getal identificeren. 
Onderstel het antwoord is 3 goed 2 raak. Het getal kan geweest zijn: 


4098 4918 4809. 


B raadt 4098 en identificeert het getal. 
Antwoord 3 goed 3 raak. Het getal was 4089 of 4819. B raadt 4089. 
Antwoord 2 goed 0 raak. Het getal kan geweest zijn 


9508 9580 8590 1598 9581 8591. 


B raadt 9581. 
Antwoord 2 goed 1 raak. Het getal kan geweest zijn 


9518 1589 8509. 


B raadt 9518. 

Antwoord 2 goed 2 raak. Het getal was 8519. 

Geval b wordt op dezelfde manier behandeld als geval a, maar levert iets minder 
mogelijkheden. Het is daarom eenvoudiger. De lezer kan het zelf zonder moeite ver- 
der uitwerken. 

Liefhebbers zouden nu kunnen nagaan, wat moet gebeuren, als de tweede infor- 
matie (bij 4567) geringer was geweest en het antwoord was geweest b.v. 1 goed 0 
raak. Ik vermoed, dat dan zes keer raden steeds uitkomst kan bieden. En zo zou 
men verder kunnen gaan, 
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